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By means of the adiabatic self-consistent collective coordinate method and the pairing- 
plus-quadrupole interaction, we have obtained the self-consistent collective path connecting 
^ ■ the oblate and prolate local minima in ®*Se and ^^Kr for the first time. The self-consistent 

' collective path is found to run approximately along the valley connecting the oblate and 

\^ , prolate local minima in the collective potential energy landscape. This result of calculation 

. clearly indicates the importance of triaxial deformation dynamics in oblate-prolate shape 

' coexistence phenomena. 

^ , §1. Introduction 
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Microscopic description of large amplitude collective motion in nuclei is a long- 
standing fundamental subject of nuclear structure physics .^''^'1^''^''^ In spite of the 
^ \ steady development in various theoretical concepts and mathematical formulations 

Ch ■ of them, application of the microscopic many-body theory to actual nuclear phe- 

nomena still remains as a challenging subject.^)"^^) Shape coexistence phenomena 
are typical examples of large amplitude collective motion in nuclei, and active inves- 
tigations are currently going on both experimentally and theoretically. ^^^"^"^^ We are 
particularly interested in the recent discovery of coexisting two rotational bands in 
^^Se and ^^Kr, which are associated with oblate and prolate intrinsic shapes.^'^ 
Clearly, these data strongly call for further development of the theory which is able 
to describe them and renew our concepts of nuclear structure. From the viewpoint 
of the microscopic mean-field theory, coexistence of different shapes implies that 
different solutions of the Hartree-Fock-Bogoliubov (HFB) equations (local minima 
in the deformation energy surface) appear in the same energy region and that the 
nucleus exhibits large amplitude collective motion connecting these different equi- 
librium points. The identities and mixings of these different shapes are determined 
by the dynamics of such collective motion. 

On the basis of the time-dependent Hartree-Fock (TDHF) theory, the selfcon- 
sistent collective coordinate (SCC) method was proposed as a microscopic the- 
ory of such large amplitude collective motions. It was extended to the case of 
time-dependent HFB (TDHFB) including the pairing correlations,!^ and has been 
successfully applied to various kinds of anharmonic vibration and high-spin rota- 
tions.l^'59i jo,),,6ii,62^i,M),.Mi,6Li,MhJli>68J,^ In order to apply this method to shape 
coexistence phenomena, however, we need to further develop the theory, since the 
known method of solving the basic equations of the SCC method, called ry-expansion 
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assumes a single local minima whereas several local minima of the po- 
tential energy surface compete in these phenomena. Some years ago, we proposed 
a new method of describing such large-amplitude collective motion, called the adia- 
batic self-consistent collective coordinate (ASCC) method.^ This method provides 
a practical scheme of solving the basic equations of the SCC methocP^l using an 
expansion in terms of the collective momentum. It does not assume a single local 
minimum, and therefore it is believed to be suitable for the description of shape 
coexistence phenomena. The ASCC method inherits the major achievement of the 
adiabatic TDHF (ATDHF) methods and, in addition, enables us to include pair- 
ing correlations self-consistently. In this method, the spurious number fluctuation 
modes are automatically decoupled from the physical modes within a selfconsistent 
framework of the TDHFB theory. This will be a great advantage when the method 
is applied to realistic nuclear problems. To examine the feasibility of the ASCC 
method, in Ref. 171]), we applied it to an exactly solvable model called the multi- 
0(4) model,E2J.EiEJC3 ^hich is a simplified version of the pairing-plus-quadrupole 
(P+Q) interaction model.EHJEHZSJ 

It was shown that this method yields a faithful 
description of tunneling motion through a barrier between prolate and oblate local 
minima in the collective potential.!^ 

In this paper, we report on our first application of the ASCC method to the 
P+Q interaction model. The major task here is to develop a practical procedure for 
solving the basic equations of the ASCC method in order to obtain the selfconsis- 
tent collective path. We investigate, as typical examples, the oblate-prolate shape 
coexistence phenomena in ^^Se and and find that the self-consistent col- 

lective paths run approximately along the valley connecting the oblate and prolate 
local minima in the collective potential energy landscape. To the best of our knowl- 
edge, this is the first time that, starting from the microscopic P+Q Hamiltonian, 
the collective paths have been fully self-consistently obtained for realistic situations, 
although a similar approach to the study of large amplitude collective motion was 
recently employed by Almehed and Walet.^''^ 

This paper is organized as follows: In §2, the basic equations of the ASCC 
method are recapitulated. In §3, we present a concrete formulation of the ASCC 
method for the case of the P+Q Hamiltonian. In §4, an algorithm to solve the basic 
equations of the ASCC method is discussed. In §5, we present results of numerical 
calculation for the oblate-prolate shape coexistence phenomena in ^'^Se and ^^Kr. 
Concluding remarks are given in §6. 

A preliminary version of this work was reported in this journal.!^ 

§2. Basic equations of the ASCC method 

In this section, we summarize the basic equations of the ASCC method.ESJ The 
basic assumption of our approach is that large-amplitude collective motions can be 
described by a set of time-dependent HFB state vectors \ 4>{q,p, ^, N)) parameterized 
by a single collective coordinate q, the collective momentum p conjugate to q, the 
particle number and the gauge angle f conjugate to A^. Then, the state vectors 
can be written in the following form: 

\(^{q,p,^,N)) = e-'^^ \^{q,p,N)) = e-'^^e'P<^^''^\(^{q)). (2-1) 

Making an expansion with respect to p and requiring that the time-dependent vari- 
ational principle be fulfilled up to the second order in p, we obtain the following set 
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of equations to determine \(j){q)), the infinitesimal generator Q{q), and its canonical 
conjugate P{q): 

HFB equation in the moving frame: 

d{<j>iq)\HMiq)\HQ))=0, (2-2) 

where 

HM{q)=H-X{q)N-^Q{q) (2-3) 
represents the Hamiltonian in the moving frame. 
Local Harmonic equations in the moving frame: 

5 [HM{q),Q{q)] - ]B{q)P{q) = 0, (2-4) 



6 {^{q)\ [HM{q), ]p{q)]-C{q)Q{q)-^^[[HM{q), {H-X{q)N)A],Q{q)]-^N 

(2-5) 

where 

B{q) = - {cPiq)\ [[H, Q{q)],Qiq)] \Hq)) ■ (2-6) 
represents the inverse mass, 

X _ 1 dBdV 

^^'^''W^Wq)^^ ^ ' 

the local stiffness, and {H — XN)a denotes the two-quasiparticle creation and anni- 
hilation parts of {H — XN) . 

The infinitesimal generators, Q{q) and P{q), satisfy 

Canonical variable condition : 

{Hq)\[Q{q),P{q)mq)) = i- (2-8) 

Once \<l){q)) and the infinitesimal generators are determined for every values of q, 
we obtain the collective Hamiltonian 'H{q,p) = \B{q)p^ + V{q) with the collective 

potential F(g) = ((^((?)|i?|<A(g)). 

In the above equations, no distinction is made between protons and neutrons 
for simplicity of notation. In the actual calculations described below, however, we 
shall explicitly treat both the neutron number A'^ and the proton number Z. 

§3. Application of the ASCC method to the P+Q model 

3.1. The P+Q Hamiltonian and signature quantum number 

Let us start with the well-known P-|-Q Hamiltonian,'^^'!^'^^ 

2 

H = Y. - E ^ [^rAr + A^Al) - I ^ dIj,D,k. (3-1) 
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^2i^= E E^25(^o4q (3-2) 



T=n,p kler 

Here, £>^^(A:/) = (A;| r^l2_ft: 10; ^'^ and x denote the pairing and quadrupole force 
strengths, respectively; c|. and are nucleon creation and annihilation operators in 
the single-particle state k, while c[ and cr denote those in the time reversed state 

k ^ 

of k. The index r indicates protons (r = p) and neutrons(r = n). Although we 
shall not explicitly mention below, it should be keep in mind that the single-particle 
index k actually includes the index r. The notation S' in the pair operators, 
and At, indicates the sum over the pairs {k,k). The factors, = {2Z/Afl'^ and 
Up = (2N / A)'^/'^ , attached to the quadrupole matrix elements guarantee equivalent 
root-mean-square radii for protons and neutrons. Following Baranger and Kumar,^ 
we take into account two major shells as a model space, and multiply a reduction 
factor Q = {Nl + 3/2) /{Nl + 5/2) to the quadrupole matrix elements D^^{kl) of 
the upper harmonic-oscillator shell, being the total oscillator quanta of the lower 
shell. According to the conventional prescription of the P-l-Q interaction,^'^'^^ 
we ignore the exchange (Fock) terms. Namely, we employ the Hartree-Bogoliubov 
(HE) approximation throughout this paper. 
We introduce the notations. 



2 

= {4±),4±),Z)£\z)(^\l)g)} {s = l- 5), (3-3) 
and write the P-l-Q Hamiltonian in the following form: 

H = Y. ^^4c. - E ^n^^F^^) + Y: ^f(-^f(-\ (3-4) 

k s=l s=l 

where Ks = {2Gn,2Gp,x,2Xj2x} for s = 1 — 5. Our Hamiltonian is invariant 
with respect to the rotation vr about the x axis. The symmetry quantum number 
associated with it is called signature, r = e~™°'. To exploit the signature symmetry, 
it is convenient to use nucleon operators with definite signatures defined by 



dk = —{ck + cj^), r = -i (q = 1/2), 

dk = ^(Cfc - Cfe), r = +i {a = -1/2), (3-5) 



1 

3^ 
1 

and their Hermite conjugates, and d^. The operators Fg^'^ are then classified 
according to the signature quantum numbers, r = ±1 (a = 0, 1): 

{^ri^\ F>2o \ F>2i \ ^22 ^} ~ +!)> 

{Dit\D!£^} (r = -l). (3-6) 

Note that D^^o^ = 0. The HB local minima corresponding to the oblate and prolate 
equilibrium shapes possess positive signature, r = +l{a = 0). Therefore, the oper- 
ators, Q{q) and P{q), generating large amplitude collective motions associated with 
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these shapes also possess positive signature. In other words, the negative signature 
degrees of freedom are exactly decoupled from the large amplitude collective motion 
of interest, so that we can ignore them. Also, it is readily confirmed that the K = 1 
components associated with the quadrupole operator ^ exactly decouple from 
the K = and 2 components in the local harmonic equations, (2.2) and (2.4). As is 
well known, they are associated with the collective rotational motions, and the large 
amplitude shape vibrational motions under consideration are exactly decoupled from 
them in the present framework. We note, however, that it is possible, by a rather 
straightforward extension, to formulate the ASCC method in a rotating frame of 
reference. By means of such an extension, we shall be able to take into account the 
coupling effects between the two kinds of large amplitude collective motion. It is cer- 
tainly a very interesting subject to study how the properties of the large-amplitude 
shape vibrational motions change as a function of angular momentum, but it goes 
beyond the scope of this paper. We note, however, that some attempts toward this 
subject was recently made by Almehed and Walet.l^ 

Thus, only the components, D^^}, are pertinent to the shape 

coexistence dynamics of our interest. They all belong to the positive signature sector, 
and we can adopt a phase convention with which single-particle matrix elements of 
them are real. In the following, we assume that this is the case. 

3.2. Quasiparticle-random-phase approximation (QRPA) at the HB local minima 

As discussed in the introduction, shape coexistence phenomena imply the exis- 
tence of several local minima in the deformation energy surface, which are solutions 
of the HB equations. Let us choose one of them and denote it as |i;^>o). The HB 
equation is written 

5{(t)o\H-Y,KNM=^, (3-7) 

T 

where A,- denote chemical potentials for protons(r = p) and neutrons(r = n). The 
quasiparticle creation and annihilation operators, ajt and a^, associated with the HB 
local minimum are defined by a^|</>o) = 0. Similar equations hold for their signature 
partners fl. They are introduced through the Bogoliubov transformations: 



k 



and their Hermite conjugate equations. (Here and hereafter, we do not mix protons 
and neutrons in these transformations.) In terms of two quasiparticle creation and 
annihilation operators, 

the RPA normal coordinates and momenta describing small amplitude vibrations 
about the HB local minimum |0o) are written 

Qp = Y,QU'^U + A^,), (3-10) 
^p = ^E^m-(4p-^mp)> (3-11) 



where the sum is taken over the proton and neutron quasiparticle pairs {p.i'), and 
p labels the QRPA modes. The amplitudes, Q'^u and P^^, are determined by the 
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QRPA equations of motion, 

6{M[H-Y^ XrNr, Qp] - -BpPp |(/.o) = 0, 



<5 (</>0| [i? - ^ XrNr, ]Pp] - CpQp |</.o> = 0, 



(3-12) 
(3-13) 



and the orthonormalization condition, (0o| [QpiPpi] — I'Up^p'. 

3.3. The HB equation and the quasiparticles in the moving frame 

For the P+Q Hamiltonian, the HB equation (|2-2j) determining the state vector 
\(p{q)) away from the local minimum reduces to 

HHQ)\hM{q)\(t>{Q))=0, 
where hM{q) is the mean-field Hamiltonian in the moving frame, 



dV_ 
dq 



h{q) = ^kidUk + did-,) - ^ {ct>{q)\ 



The state vector \(f){q)) can be written as a unitary transform of 



(3-14) 

(3-15) 
(3-16) 



p,i/ 1 ) 



(3-17) 



pu 



where the sum is taken over the proton and neutron quasiparticle pairs {fJ-i'). The 
quasiparticle creation and annihilation operators, ap{q) and ap{q), associated with 
the state \(p{q)), which satisfy the condition, ap{q)\(p{q)) = 0, are written: 



(3-18) 



at (g) EE e''((?)at e-^(^) = J] (c/^.((?)4 + Vp,{q)a 

V 

a-p{q) = e\q)a-pe-'^'i^ = (v-p,{q)o^l + U-p,{q)aA, 



where 



Upv(q) Vpv{q) 

Vpuiq) Upp{q) 



V 



sin(\/^) \ 



sm 



(3-19) 



Here, 9 in the r.h.s. denotes the matrix composed of Opuiq), and it is understood 
that its elements corresponding to those in the l.h.s. should be taken. 

In terms of the quasiparticle operators defined above, the mean-field Hamilto- 
nian in the moving frame hmiq), the neutron and proton number operators, Nr, 
and the pairing and quadrupole operators, Fg^^ , are written in the following form: 



/^M(g) = (</'(g)|/iM(g)|</'('7)) + (3-20) 
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(3-21) 



= {cP{q)\ Fi±) 10(g)) +^F^^\^,u)[Al,{q) ± A,,{q)) 



+ J2^its(H{Bf.u{q) + B^,iq)), 



(3-22) 



where 



Alp{q) ^ al{q)al{q), A^p{q) = ap{q)a^{q), B^^{q) = al{q)a„{q). (3-23) 

Note that Ep_{q) = E^{q) and also that the equahties, Fj^^{p,r') = F^l{^v), hold for 
the operators under consideration. Explicit expressions for the expectation values 
and the quasiparticle matrix elements appearing in the above equations are given in 
Appendix A. 

3.4. Local harmonic equations in the moving frame 

We can represent the infinitesimal generators, Q{q) and P{q), in terms of Aj)p(g) 
and A^p{q) as 



where the sum is taken over the proton and neutron quasiparticle pairs (^z^). For 
the P+Q Hamiltonian, the local harmonic equations, ()2-4() and ()2-5|) . in the moving 
frame reduce to 



HHq)\ [hM{Q),QiQ)] - E/fc^i"^ - -B{q)P{q) IHq)) = 0, (3-26) 



S {m\ [hM{q) , -B{q)P{q)] - ^ /pTi " B{q)C{q)Q{q) - ^ /^^^^ 




(3-24) 



(3-25) 



s 



s 



s 




(3-27) 



s 



T 



where the quantities /q j etc. are given by 



fQ,s ^ {Hq)\ [Fs^-\QiQ)] mi)) = 2K,{F(-\Q{q)), (3-28) 

^ (0(g)| ]Biq)P{q)] |0(g)) = 2K,i?(g)(Fi+), P(<?)), (3-29) 

/ff = -l^s {m\ [Ri^\Q{q)] \m) = Ks{Ri+\Qiq)), (3-30) 

f^^,^B{q)-^. (3-31) 
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Here we have introduced the notations 



where {h{q) — Y^^ Xr{q)NT)A represents the A^^p{q) and A^p{q) parts of the operator 
in the parenthesis. We also use the notations, 

(FtlQiQ)) ^Y.^s'\}^^)QAl)^ etc. 



(3-33) 



Note that /q I, and are hnear functions of Qf_ip{q) or Pi_ip{q). 

One can easily derive the following expressions for the matrix elements, Qfj,u{q) 
and Pf^p{q), from the local harmonic equations in the moving frame, H3-26() and 

s s 

+ Ri-\f,i?)f^Q-l +Y.N^(i,u)fN,r} (3-34) 

r 



+ u;\q)Y^g2{^^u)F(-H^,u)fi^-l 



where fi^^ = + /^^, 



(+) 



and 



(3-35) 



(3-36) 



{E^+Epr-u;^{qy "^'"^^ ' {E ^ + Epf - Lo^qY 

Note that w^, representing the square of the frequency uj{q) = ^/B(q)C{q) of the 
local harmonic mode, is not necessarily positive. The values of B{q) and C{q) depend 
on the scale of the collective coordinate q, while Lo{q) does not. Namely, the scale 
of q can be chosen arbitrarily without affecting the frequency uj{q). We thus require 
B{q) = 1 everywhere on the collective path to uniquely determine the scale of q. 

Inserting expressions (3.34) and (3.35) for Q^p{q) and P^p{q) into Eqs. (3.28)- 
(3.30), and combining them with the orthogonality condition to the number opera- 
tors, 

{4>{q)\[Nr.P{q)]\m) = 2i{P{q),Nr) = 0, (3-37) 
we obtain linear homogeneous equations: 



PR,Q 



qQ,PR 



PR,PR 



\ 



qPR,N 



1 nN,C 

p(-) A+) 



S. 



N,PR 



S. 



N,N 



( \ 



A+) 

•I PR,! 



\ fNy J 



(3-38) 



for the vectors, /g , f p^^, and defined by 



J Q — ^Jq,i ' Jq,i J' 

A+) ^ J A+) A+) A+) A+) X 

JPR—IJPR,1^ J PR,2' JPR,V JPR,5J^ 
/at = {fN,n, fN,p}- 



(3-39) 

(3-40) 
(3-41) 
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Here, 

Sff = 2{Ft\F^/\, + 2(Fi-),iiir^),, - ^dsg,, (3-42) 



sQfR^2iF(-\F^/\„ (3-43) 
Sff^2{F(-\N,,)g,, (3-44) 

+ iRM,F^~\,, (3-45) 
^P^,P^ . 2(fW,fW),, + {R^\F^/\, - (3.46) 

Kg 

5f,!^'^ ^ 2(Fi+),iV^05i + (M+\A^r')p2, (3-47) 

^ ^ (iV,,i?(7^),,, (3-48) 

5f/^^(iV.,Fi,+^),„ (3-49) 

5f/^(iV„iV^,)9i, (3-50) 
with the notations, 

iF^~\F^-^)g,^^F^-\^i?)g,ifii?)F^~\^i?), etc. (3-51) 

Equation (|3-38)) is of the form 

8 

J2Saa'{u^HQ))U' =0, (3-52) 

(t'=1 

for the vector / composed of 



.(-) A+) 

Q ' J PW 

,{-) A~) A+) A+) A+) A+) 



{/.=i_8} = {fQ\ f^i fN} (3-53) 



— {/q,1' /q,2' /p_R,l' /pp,2' /pp,3' /pR,5' /a^,"' /a^.p}' (3-54) 

SO that the frequency co of the local harmonic mode is determined by the condition 
det S = 0. The normalizations of / are fixed by 

{4>{q)\[Q{q),P{q)Mq)) = 2i{Q{q), P{q)) = i. (3-55) 

Note that lo^ represents the curvature of the collective potential, 

(3-56) 

for the choice of coordinate scale with which B{q) = 1. 

In concluding this section, we mention that the reduction of the local harmonic 
equations to linear homogeneous equations like (|3-52() can be done for any effective 
interaction that can be written as a sum of separable terms. Below, we call the local 
harmonic equations in the moving frame "moving frame QRPA" for brevity. 

§4. Procedure of calculation 

4.1. Algorithm to find collective paths 

In order to find the collective path connecting the oblate and prolate local min- 
ima, we have to determine the state vectors \(t>{q)) and the infinitesimal generators, 
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Q{q) and P{q), by solving the moving frame HB equation (|3-14|) and the moving 
frame QRPA equations, (|3-26|) and (|3-27|) . Since Q{q) depends on \(j){q)) and vice 
versa, we have to resort to some iterative procedure. We carry out this task through 
the following algorithm. 

Let us assume that the state vector, \(j)iq)), and the infinitesimal generators, Q{q) 
and P{q), are known at a specific point of q. We then find the state vector, \4>{q+6q)), 
and the infinitesimal generators, Q{q + 6q) and P{q + 5q), in the neighboring point, 
q + 5q, through the following steps. 

Step 1: Construct a state vector at the neighboring point, q + 5q, with the use of 

+ 5(z))(°) = e-^^'?^(^)|,/.((?)). (4-1) 

Though \(j){q + 5q))^^^ does not necessarily satisfy the moving frame HB equation, 
(|3-14|) . we can use this state vector as an initial guess at g + 5q. 
Step 2: Solve the moving frame HB equation ()3-14() with the use of Q(0)(g + (5(7) = 
Q{q) as an initial guess for Q{q + 5q), and obtain an improved state vector \(l){q + 
6q))^^\ In doing this, we find it important to impose the constraint, 

{4>{q + Sq)\Q{q)\4>{q + dq)) = 6q (4-2) 

for the increment 6q of the collective coordinate q, together with the constraints, 

{^{q + 6q)\Nr\HQ + SQ)) =Nr, r=p,n (4-3) 

for the proton and neutron numbers {Np = Z, Nn = N). The constraint ()4-2() is 
easily derived by combining Eq. (|4-lj) with the canonical variable condition (|2-8p . 
Details of this step is described in Appendix B. 

Step 3: Solve the moving frame QRPA equations, H3-26() and H3-27|) . with the use of 

\(j){q + (5g))(^) to obtain Q'^^\q + 5q) and P'^^\q + 6q). 

Step 4: Go to Step 2 and solve Eq. ^^T^i with the use of Q^^\q + 5q). 

If the iterative procedure. Steps 2-4, converges, we obtain the selfconsistent 
solutions, Q{q + 5q), P{q + 5q) and \(l){q + 5q)), that satisfy Eqs. ()3-14|) . H3-26() and 
H3-27() simultaneously at q + 5q. Then, we go to Step 1 to construct an initial guess 
\(f){q + 26q))^^^ for the next point, q + 26q, and repeat the above procedure. In this 
way, we proceed step by step along the collective path. 

The above is a brief summary of the basic algorithm. In actual numerical 
calculations, we start the procedure from one of the HB local minima and choose the 
lowest frequency QRPA mode as an initial condition for the infinitesimal generators, 
Q and P, at q = 0. Under ordinary conditions, we can proceed along the collective 
path following the procedure described above. In some special situations, however, 
we need additional considerations concerning the choice of the initial guess, Q^^\q-\- 
Sq), in Step 2. Actually, we meet such situations in some special regions of the 
collective path in ^■^Kr. We shall give a detailed discussion on this point in §§ 5.3. 

We have checked that the same collective path is obtained by starting from the 
other local minimum and proceeding in an inverse way. 

4.2. Details of calculation 

In the numerical calculation, we use the spherical single-particle energies of 
the modified oscillator model of Ref. I83j) . which are listed in Tabled and follow 
the conventional prescriptions of the P+Q interaction model,!^ except that the 
pairing and quadrupole interaction strengths, Gr and are chosen to approximately 
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reproduce the pairing gaps and quadrupole deformations obtained in the Skyrme- 
HFB calculation by Yamagami et o/P* Their values are: G„ = 0.320 (0.299) 
Gp = 0.320 (0.309) and x' = X^^ = 0.248 (0.255) in units of MeV for ^^Se f ^Kr), 
where b is the length parameter given by 62 ^ 4^)1/3^2^1/3^ ^he pairing gaps, 
A^-^p^n, and deformation parameters, (3 and 7, are defined as usual through the 
expectation values of the pairing and quadrupole operators: 

Z\.(g)=G.(0(g)|j;44|</.(g)), 

fcGr 

Pcosj = x'{HQ)\DiV\HQ))/{f^ob^), 
Psinj = V2x'{HQ)\DiP\HQ))/{t^ob^), 
where hujo denotes the frequency of the harmonic oscillator potential. 

Table I. Spherical single-particle orbits and their energies used in the calculation. The energies 
relative to those of lgg/2 are given in units of MeV. 



orbits 


1/7/2 


2p3/2 


1/5/2 


2pi/2 


I59/2 


2^5/2 


lfl7/2 


3Sl/2 


2rf3/2 


protons 


-8.77 


-4.23 


-2.41 


-1.50 


0.0 


6.55 


5.90 


10.10 


9.83 


neutrons 


-9.02 


-4.93 


-2.66 


-2.21 


0.0 


5.27 


6.36 


8.34 


8.80 



§5. Results of calculation 

5.1. Properties of the QRPA modes at the local minima in ^^Se and '^^Kr 

For the P-l-Q Hamiltonian described in § 4, the lowest HB solution corresponds 
to the oblate shape, while the second lowest HB solution possesses the prolate shape, 
for both ^^Se and ^^Kr (see Table|ll|. Their energy differences are 0.30 and 0.82 MeV 
for ®®Se and ''^Kr, respectively. In the QRPA calculations at these local minima, we 
obtain strongly collective quadrupole modes with low- frequencies. They correspond 
to the P and 7 vibrations in deformed nuclei with axial symmetry. Although the 
former in fact contains pairing vibrational components, we call it /3 vibration because 
the transition matrix elements for the quadrupole operator are enhanced. (A 
neutron pairing vibrational mode appears as the second QRPA mode at the oblate 
minimum in ^^Kr; see Table ??.) Let us note that there is an important difference 
between ^^Se and ^^Kr concerning the relative excitation energies of the f3 and 7 
vibrational modes: In the case of ^®Se, the frequencies of the 7 vibrational QRPA 
mode are lower than those of the /3 vibrational one for both the oblate and prolate 
local minima. The situation is opposite in the case of ^^Kr, namely, the /3 vibrations 
are lower than the 7 vibrations. As we shall see in the succeeding subsections, this 
difference leads to an important difference in the properties of the collective path 
connecting the two local minima. 

5.2. Collective path connecting the oblate and prolate minima in ^^Se 

As the 7 vibrational mode is the lowest-frequency and most collective QRPA 
mode at the prolate local minimum, we have chosen this mode as an initial condition 
for solving the basic equations of the ASCC method, and carried out the procedure 
described in §§ 4.1. We thus obtained the collective path connecting the oblate and 
prolate local minima in ^'^Se, which is plotted in Fig. 1(a). As we have extracted 
the collective path in the TDHB phase space, which has huge degrees of freedom. 



(4-4) 

(4-5) 
(4-6) 
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Table II. The equilibrium quadrupole deformation parameters (/3, 7), the pairing gaps (At) in 
units of MeV, the QRPA eigenenergies ftLJp=i,2 in units of MeV, and the relevant quadrupole 
transition matrix elements squared, |A/pp = | (pli'gj^-' |0)p (p — 1,2). Here, \p) and |0) denote 
the QRPA one-phonon and the ground states. The notations, (/3), (7), and (/in), in the eighth 
column, respectively, indicate the /3-, the 7- and the neutron-pairing vibrational modes; the 
|Mpp values for K = 0, 2, and are presented, respectively, in the Weisskopf units. 





13 


7 


A„ 


Ap 




\Mi\'^ 


UJ2 


IM2I' 


'=«Se (prolate) 


0.234 


0° 


1.34 


1.42 


1.02(7) 


33.66 


1.91(/3) 


12.19 


®*Se (oblate) 


0.284 


60° 


1.17 


1.27 


1.55(7) 


13.64 


2.25(/3) 


7.67 


''^Kr (prolate) 


0.376 


0° 


1.15 


1.29 


1.60(/3) 


12.97 


1.67(7) 


14.61 


''^Kr (oblate) 


0.354 


60° 


0.86 


1.00 


1.15(/3) 


5.37 


1.91(Zi„) 


0.19 



the path drawn in this figure should be regarded as a projection of the cohective 
path onto the (/3,7) plane. Roughly speaking, the collective path goes through the 
valley that exists in the 7 direction and connects the oblate and prolate minima. 
If f3 is treated as a collective coordinate and the oblate and prolate shapes are 
connected through the spherical point, the variation of the potential energy would 
be much greater than that along the collective path we obtained. The potential 
energy curve V{q) along the collective path evaluated using the ASCC method is 
shown in Fig. 1(b). Since we have defined the scale of the collective coordinate q 
such that the collective mass is given by M{q) = B{q)-^ = 1 MeV'^, the collective 
mass as a function of the geometrical length s along the collective path in the (/?, 7) 
plane can be defined by 

M(.(g)) = M(g)(^)-2 (5-1) 

with ds'^ = dp'^ + /3^d7^. This quantity is presented in Fig. 1(c) as a function of 
q. The triaxial deformation parameter 7 is plotted as a function of q in Fig. 1(d). 
Variations of the pairing gaps, A^-iq), and of the eigen- frequencies of the moving 
frame QRPA equations along the collective path are plotted in Figs. 1(e) and 1(f). 
The solid curve in Fig. 1(f) represents the frequency squared, u)'^{q) = B{q)C{q), 
given by the product of the inverse mass B{q) and the local stiffness C{q), of the 
solutions of the moving frame QRPA equations, that corresponds to the 7- vibration 
in the oblate and prolate limits. These QRPA solutions determine the infinitesimal 
generators Q{q) and P{q) along the collective path. For reference sake, we also 
present in Fig. 1(f) another solution of the moving frame QRPA equations, which 
possesses the /3-vibrational character and is irrelevant to the collective path in the 
case of ^^Se. Note that the frequency of the 7-vibrational mode becomes imaginary 
in the region, 12° < 7 < 45°. These results should reveal interesting dynamical 
properties of the shape coexistence phenomena in ^^Se. For instance, the large 
collective mass in the vicinity of 7 = 60° (Fig. 1(c)) might increase the stability of 
the oblate shape in the ground state. Detailed investigation of these quantities as 
well as solutions of the collective Schrodinger equation will be given in the succeeding 
paper.l^ 

5.3. Collective path connecting the oblate and prolate minima in '^'^Kr 

In contrast to ®^Se, the lowest-frequency QRPA mode is the /3 vibration at 
the prolate local minimum in ^^Kr. Therefore, we have chosen this mode as the 
initial condition at the prolate minimum and started the procedure of extracting the 
collective path. Then, the collective path first goes in the direction of the /? axis on 
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the (/9, 7) plane. As we go along the /? axis, we eventually encounter a situation that 
the two solutions of the moving frame QRPA equations compete in energy, and they 
eventually cross each other. Namely, the character of the solution with the lowest 
value of (j^ = BC changes from the (3 vibrational to the 7 vibrational ones at some 
point on the collective path. If only the solution Qi (q) with the lowest value of at 
the previous point q is always chosen as an initial guess for Q{q + Sq) in Step 2 of the 
algorithm described in §§ 4.1, then the direction of the collective path on the {P,j) 
plane change abruptly from the /3 to the 7 directions immediately after the crossing 
point (in the vicinity of the point C in Fig. 3 presented below), and the numerical 
algorithm outlined in §§ 4.1 fails at this point: During the iterative procedure of 
solving the moving frame HB equation, we encounter a situation where the overlap 
(Qiq), Qiq + Sq)) between the infinitesimal generators Q at the neighboring points, 
q and q + 6q, vanishes, because the former possesses K=0 whereas the latter has 
K=2. The numerical algorithm (details of which is described in Appendix B) then 
stops to work just at this point, where the overlaps (Nr, Q{q + Sq)) also vanish by 
the same reason. This problem occurs even if wc decrease the step size 5q. We 
find, however, that we can avoid this difficulty by employing a more suitable initial 
guess for Q{q + 5q). Namely, we take a linear combination of the two solutions, 
Qi{q) and Q2iq) at the previous point q, Q^°\q + 5q) = (1 - e)Qi{q) + eQ2{q), 
with a small coefficient £, as an initial guess. This improvement is just for starting 
the iterative procedure at the next point, q + 5q, on the collective path, so that 
the self-consistent solution, Q{q + 5q), obtained upon convergence of the iterative 
procedure, of course, do not depend on the values of e. For instance, wc obtain an 
axially symmetric solution \(p{q + 6q)) and a generator Q{q + 6q) preserving the K 
quantum number in the region with /3 > 0.24 around the prolate minimum, even 
when we start the iterative procedure using an initial guess for Q{q + 5q) that breaks 
the axial symmetry. We confirmed that this is indeed the C3jSG ctS long SjS £ is 3i small 
finite value around 0.1. This special care is needed only near such crossing points (as 
shown below in Figs. 3-5) where two solutions of the moving frame QRPA equations 
with different K quantum numbers compete in energy. 

With the improved algorithm mentioned above, we have successfully obtained 
the smooth deviation of the direction of the collective path from the (3 axis toward 
the 7 direction (see Fig. 2(a)). We note that the character of the lowest uP solution 
of the local harmonic equations also gradually changes from the (3 vibrational to the 
7 vibrational ones (see Fig. 2(f)). The details of the turn over region is presented in 
Fig. 3. One can clearly see in Fig. 3(a) a gradual onset of axial-symmetry breaking in 
the solutions \4>{q)) of the moving frame HB equation. One can also see in Fig. 3(b) 
an avoided crossing between the lowest two solutions of the moving frame QRPA 
equations associated with mixing of the components with K={) and 2. After the 
smooth turn to the 7 direction, the 7 value increases keeping the /3 value roughly 
constant, and the collective path eventually approaches to the 7 = 60° axis. Then, 
we again encounter a similar situation. Adopting the improved algorithm, wc have 
confirmed that the character of the lowest solution smoothly changes, this time, 
from the 7 vibrational to the (3 vibrational ones. The collective path thus merges 
with the 7 = 60° axis, and finally reaches the oblate minimum. 

We have also carried out the calculation starting from the oblate minimum and 
proceeded in an inverse way, and obtained the same collective path. This should 
be regarded as a crucial test of consistency of our calculation. Figure 4 shows the 
details of this test: The collective path that started from the prolate minimum and 
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turned into the 7 direction gradually merges with the 7 = 60° axis. On the other 
hand, moving in the opposite direction, we see a gradual onset of axial-symmetry 
breaking in the collective path that started from the oblate minimum. We see that 
the two results of calculation for the collective path nicely agree with each other. 
The importance of taking account of the mixing between the (3- and 7-vibrational 
degrees of freedom in solving the moving frame HB and QRPA equations is again 
demonstrated in Fig. 5, which shows the details of the turn over region from the 
7 = 60° axis. 

Although the collective path drawn in Fig. 2(a) should be regarded as a projec- 
tion of it onto the (/?, 7) plane, the result of calculation indicates that the collective 
path runs roughly along the valley in this plane. The potential energy curve V{q), 
the collective mass M(s{q)), variations of the pairing gaps, Ar^q), are presented in 
Figs. 2(b), (c), (e), respectively. Their properties are similar to those for ^*Se. In 
particular, we notice again a significant increase of M(s(g)) in the vicinity of the 
oblate minimum. 

Quite recently, Almehed and Walet studied the oblate-prolate shape coexistence 
phenomenon in ^^Kr by means of an approach similar to the ASCC method but with 
some additional approximations,!^ and found a collective path going from the oblate 
minimum over a spherical energy maximum into the prolate secondary minimum. 
We have also obtained such a collective path when we impose the axially symmetry 
to the solutions \4>{q)) of the moving frame HB equation and always use only K=0 
solutions of the moving frame QRPA equations. But, when we relax such symmetry 
restrictions and follow the lowest u)"^ solution of the moving frame QRPA equations, 
we obtain the collective path presented in Fig. 2, which breaks the axial symmetry. 
The reason of this disagreement is not clear at present. With their parameters of 
the P+Q Hamiltonian, they did not encounter the character change of the lowest 
moving-frame-QRPA mode on the collective path, from the P vibrational to the 7 
vibrational ones. However, it is interesting to observe that they in fact encountered 
the avoided crossing with a 7-vibrational mode, similar to the one shown in Fig. 2(f), 
and obtained a collective path that turns into the triaxial plane in their calculation 
for states with angular momentum 1 = 2. 

§6. Concluding Remarks 

We have applied the ASCC method to the oblate-prolate shape coexistence phe- 
nomena m •^^Se and ^^Kr. It was found that the self-consistent collective paths run 
approximately along the valley connecting the oblate and prolate local minima in the 
collective potential energy landscape. This is the first time that the self-consistent 
collective paths between the oblate and prolate minima have been obtained for re- 
alistic situations starting from the microscopic P-|-Q Hamiltonian. Recently, the 
generator coordinate method has often been used to describe a variety of shape 
coexistence phenomena, with /? employed as the generator coordinate.'^''^ The 
triaxial shape vibrational degrees of freedom are also ignored in the extensive vari- 
ational calculations by Tiibingen group.'^''^ The result of the ASCC calculation, 
however, strongly indicates the necessity of taking into account the 7 degree of free- 
dom, at least for the purpose of describing the oblate-prolate shape coexistence in 
^®Se and ^■^Kr. In order to evaluate the mixing effects between the oblate and prolate 
shapes taking into account the triaxial deformation dynamics, we have to quantize 
the classical collective Hamiltonian obtained in this paper and solve the resulting 
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collective Schrodinger equation. This will be the subject of the succeeding paperP^ 
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Appendix A 

Explicit expressions of the quasiparticle matrix elements 

Combining the successive Bogoliubov transformations, (3.8) and (3.18), the 
quasiparticles, ali{q) and af^{q), associated with the state |(/'((?)), can be written 
in terms of the nucleon operator, d\. and d^, as 

Making use of the inverse transformation, 

\d-k ) ^\ V-,^{q) U-,^{q) 

one can easily derive the exphcit expressions of the expectation values and the matrix 
elements of the operators, Fs^\ appearing in Eq. (|3-22j) : 

{Hq)\F^s=i2\Hq)) = -2 j;5;(A;fe|4t)„_p|0)t/,^(g)n^(g), 

M k 

k 

Fb}=iM) = Y.^kk\A^^X,p\Q){UkMy-ku{Q) ± UUlWi.iq)), 

k 

m)\ m)) = '^Y.Y.^^\^iK=o,2\m-Myi-M^ 

kl 

kl 

4i3,5(^'^) = Y.^k\D^A=o,2\l){ykMUUq) - VuMVi-^iq)), (A-3) 
kl 

The expectation values of the anti-Hermitian operators, ^r=n,p, vanish. The quan- 
tities, {k\D^\l), etc., appearing in the above expressions denote matrix elements 
between the single-particle states defined by Eq. (|3-5|) : 

(fc|D(+)|/)^(o|4i)^+)dl|o), 

(fcfc|4±)|0) ^ (01444^^10), etc., (A-4) 

where |0) is the vacuum for the nucleon operators {d'^^d). The matrix elements of 
the Bogoliubov transformations, (A.l) and (A. 2), possess the following symmetries: 

U,k = U-^k = Uk, = U-,f„ V^-, = -V-,k = Vk-, = -V-u^. (A-5) 
It is also easily seen that equalities, 

(^l^2S=0,2lO = (fc|^S=0,2lO, (A-6) 

fW(/2P) = f£)(/xz.), F'i'}{v^^) = ±F'i}{^lv), (A7) 
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hold for the pairing and quadrupole operators, {A^\a^\d^\d^^^}, under con- 
sideration. The expectation values, ((^(g)] Nt \4'(,q))j a-nd the matrix elements, A''^(/u) 
and NB^rifJ-), of the neutron and proton number operators are readily obtained from 
those of F^t^l by replacing {k\D^^\l) with 6ki, and restricting the sum over the 
single-particle index k to neutrons or protons. 

Appendix B 

Solving the moving frame HB equation 



We solve the moving frame HB equation using a method similar to the imaginary 
time method.!^ Let \(t)^'^\q)) be the state vector in the iterative step i. We first 
calculate the mean- field Hamiltonian associate with it: 



{q) = ^ 6fc(44 + d\d-,) - 

k s 



Using the quasiparticle operators h^i^^ and 6^*' defined by 



bf 0«(g)) = O, 

we then generate a state vector in the {i + 1) step as 

<^(i+i)(g)^ = exp 1(^+1) 

'm{q)-Y^Xi:^+^){q)N, 



(B-1) 



(B-2) 



+ E[U'\q) - X^^^'\q)Nr - fi^'+'\q)Q{q) 



X 



where e is a small parameter, 



(B-3) 



(B-4) 



and the subscripts, + and — , denote the two-quasiparticle creation and annihilation 
parts of the operator in the parenthesis, respectively. It should be noted that, in con- 
trast to the conventinal imaginary time method, the unitary operator, exp , IS 
used here so that the normalization is preserved during the iteration. The Lagrange 
multipliers \T~^^\q) and fJ.^^~^^\q) are determined by the constraint equations. 



0(^+1) (g) Qiq-6q) cl^^'+'Hq)) = 6q, 



(B-5) 



where and N^°^ are the neutron and proton numbers of the nucleus under 

consideration. Similar but slightly different constraints are utilized by Almehed 
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and Walet!^ Expanding the left-hand sides up to first order in x^^'^^^ we obtain 
equations determining them: 

iNp,Nn) 

{Q{q-5q),Nn) 



where the quantities, {N^-, N^-i), {NT-,Q{q)) etc., are defined by (3.33), except that 
the coefficients, Nr{fJ-), Qij.u{q) etc., involved in these quantities are here defined with 
respect to the two-quasiparticle creation and annihilation operators, and -X^^p. 

Using the state vector \4>^''~^^\q)) , we calculate the mean- field Hamiltonian h^''~^^\q) 
in the (i + 1) step, and repeat the above procedure until convergence is attained. 
The mean-field Hamiltonian thus obtained takes the following form: 

hM{q) = h{q) - A,(g)iV, - fi{q)Qiq) 

T 

= {ct>{q)\ hliq) \m) + E {bl{l)K{q) + hl{q)hM) ■ (B-7) 

Finally we introduce the quasiparticle operators, a^^{q) and ai_i{q), that diagonalize 
hniq)- 

huiq) = (</'(9)| /im(9) 10(g)) +J^^^(g)(aj,(<7)a^(g)+at(g)a^(g)). (B-8) 

It is easy to see that ^{q) = dV/dq. In actual calculations, the above procedure is a 
part of the double iterative algorithm described in section 4. Namely, we carry out 
the above iterative procedure using the constraint operator Q{q)^"'^ that is obtained 
in the n-th iteration step determining the infinitesimal generators, Q{p) and P{q)- 



{Q{q-6q),Np) 
( (M°^ 

(iv(°) 

{Q{q- 



iNn,Q{q)) 
{Np,Q{q)) 
{Q{q-dq),Q{q)) 



I xi:^'\q) 

(q) 

V i^^'^'Kq) 



p 



V i^q 



-{Nn)^'^)/2e + {h('\q),Nn) 
-(iV,)«)/2e + (/i«(g),iVp) 
6q))(^))/2e + {h('){q),Qiq-6q)), I 



(P-6) 
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Fig. 1. Results of calculation for ^*Se. (a) The bold curve represents the ASCC path projected 
onto the (/3, 7) plane, which connects the oblate and the prolate minima designated by filled 
circles. The contour lines are calculated using the conventional constrained HB method and 
plotted for every 50 keV. (b) Collective potential V{q) plotted as a function of the collective 
coordinate q. Here the origin of q is chosen to coincide with the prolate local minimum, and 
its scale is defined such that the collective mass is given by M{q) = 1. (c) Collective mass 
M{s{q)) with respect to the geometrical length s{q) along the collective path in the (/3, 7) 
plane, plotted as a function of q. (d) The triaxiality parameter 7 as a function of q. (e) Neutron 
and proton pairing gaps, A„ and Ap, as functions of q. (f) Lowest two eigen-frequencies squared 
(i.e., uj^ = BC) of the moving frame RPA, plotted as functions of q. These modes at triaxial 
deformed shapes are more general than the ordinary /3- and 7-vibrations at the oblate and the 
prolate limits, and contain both components. The symbols, /? and 7, are used, however, in order 
to indicate the ma.inr comnonents of the movincr frame RPA modes. 
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Fig. 2. Results of calculation for Kr. Notations are the same as in Fig. 1 except for the foUowings: 

In (a), the contour lines arc plotted for every 100 keV. In (f), the lowest three eigen-frequencies 
squared (i.e., = BC) of the moving frame RPA, are plotted as functions of q. As mentioned in 
the caption to Fig. 1, these modes at triaxial deformed shapes are more general than the ordinary 
/?- and 7-vibrations at the oblate and the prolate limits, and contain both components. The 
symbols, and 7, are used, however, in order to indicate the major components of the moving 
frame RPA modes. Likewise, the symbols An is used to indicate that the major component is 
the neutron pairing vibrational mode. 
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Fig. 3. Enlargement of the turn over region of Figs. 2(a) and 2(f) for ^^Kr, where the collective 
path turns into the 7 direction. The step size 6q = 0.0157 and e = 0.1 are used in the numerical 
calculation. Every steps in 5q arc represented by filled circles and connected by solid linos. The 
points designated A, B, C, D on the collective path in (a) correspond to those in (b) which 
displays the squared frequencies, w^, of the lowest two solutions of the moving frame QRPA 
equations as functions of the collective coordinate q. The open circles represent those obtained 
in the calculation with e = where the mixing efi'ects between the K=0 and 2 components 
are totally ignored. The points designated A', B', C in the latter calculation correspond to 
the points A, B, C in the former calculation. In the latter calculation, we cannot get the point 
corresponding to D because the problem discussed in the text occurs in the numerical algorithm. 
It was chocked that the same collective path is obtained with use of Sq = 0.0314 except that 
the distances between the successive points are doubled. 
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Fig. 4. Enlargement of the turn over region of Figs. 2(a) for ^^Kr, where the collective path (solid 
curve) coming from the prolate minimum merges with the 7 = 60° axis (dotted line). Every 
steps in 6q are represented by open circles and connected by solid lines. For comparison, the 
result of calculation starting from the oblate minimum and moving in the opposite direction is 
shown by open squares. Slight deviations from the solid curve indicate the degree of precision 
of the present numerical calculation. The step size Sq = 0.0157 and e = 0.1 are used in both 
cases. The collective path obtained by these different calculations well agree with each other. 
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Fig. 5. Enlargement of the turn over region of Figs. 2(a) and 2(f) for ^^Kr, where the collective 
path (solid curve) coming from the oblate minimum starts to deviate from the 7 = 60° axis 
(dotted line) . The numerical calculation was done starting from the oblate minimum and using 
the step size 6q = 0.0157 and e = 0.1. Every steps in 5q are represented by filled circles 
and connected by solid lines. The points designated A, B, C, D on the collective path in 
(a) correspond to those in (b) which displays the squared frequencies, w^, of the lowest two 
solutions of the moving frame QRPA equations as functions of the collective coordinate q. Note 
that the values of q in this figure are measured from the oblate minimum. The open circles 
represent those obtained in the calculation with £ = where the mixing effects between the 
K=() and 2 components are totally ignored. The points designated A', B', C in the latter 
calculation correspond to the points A, B, C in the former calculation. In the latter calculation, 
we cannot get the point corresponding to D because the problem discussed in the text occurs 
in the numerical algorithm. Slight wiggles along the successive points seen in (b) are due to 
numerical error, and they indicate the degree of precision of the present numerical calculation. 
It was checked that the same collective path is obtained with use of Sq = 0.0314 except that 
the distances between the successive points are doubled. 



